SPARTA (26 Feb 2021)
################################################################################
# beam of particles striking the surface at an inclined angle
# free molecular flow (no collisions)
#
# Note:
#  - The "comm/sort” option to the “global” command is used to match MPI runs.
# The "comm/sort" option should not be used for production runs.
################################################################################

seed	    	    123456
dimension   	    3
global              gridcut 0.0 comm/sort yes

boundary	    	oo oo oo


create_box          -11 11 -11 11 0 10
Created orthogonal box = (-11 -11 0) to (11 11 10)
create_grid 	    2 2 2
WARNING: Could not acquire nearby ghost cells b/c grid partition is not clumped (/ascldap/users/stamoor/sparta_stanmoore1/src/grid.cpp:410)
Created 8 child grid cells
  CPU time = 0.00168109 secs
  create/ghost percent = 90.668 9.33201
balance_grid        rcb cell
Balance grid migrated 4 cells
  CPU time = 0.000834703 secs
  reassign/sort/migrate/ghost percent = 71.6367 1.37104 13.9674 13.0249

global		    	nrho 1e10 fnum 1e6

species		    	air.species O CO CO2 O2 C
mixture		    	air O O2 vstream 0 1000 -1000

mixture             air O   frac 1.0
mixture             air CO  frac 0.0
mixture             air CO2 frac 0.0
mixture             air C   frac 0.0
mixture 			air O2 	frac 0.0


surf_collide        1 cll 300.0 0.5 0.5 0.5 0.5

read_surf			base_plate.surf
  8 points
  12 triangles
  -11 11 xlo xhi
  -11 11 ylo yhi
  0 1 zlo zhi
  1 min triangle edge length
  11 min triangle area
  0 0 = number of pushed cells
  4 0 = cells overlapping surfs, overlap cells with unmarked corner pts
  4 0 4 = cells outside/inside/overlapping surfs
  4 = surf cells with 1,2,etc splits
  4356 4356 = cell-wise and global flow volume
  CPU time = 0.00130057 secs
  read/check/sort/surf2grid/ghost/inout/particle percent = 16.242 12.1907 0.843263 61.6865 9.03758 15.5454 0.60495
  surf2grid time = 0.000802279 secs
  map/comm1/comm2/comm3/comm4/split percent = 24.6657 6.30015 4.10104 3.0312 6.00297 37.7117

##################################### SURF REACT ADSORB ######################################
##################################### SURF OPTION ############################################

#surf_react        	 adsorb_test_ps1 adsorb ps sample-PS_1.surf nsync 1 surf 1000 6.022e18 O CO
#surf_modify 		 all collide 1 react adsorb_test_ps1

surf_react        	adsorb_test_ps2 adsorb ps sample-PS_2.surf nsync 1 surf 1000 6.022e18 O CO
surf_modify 		all collide 1 react adsorb_test_ps2

########################## BEAM ############################################################
# Beam at multiple points so that different processors handle the surface collisions

region              circle2 cylinder z  6 -10 1 -INF INF
region              circle3 cylinder z -6 -10 1 -INF INF

fix                 in2 emit/face/file air zhi data.beam beam_area_2 nevery 100 region circle2
fix                 in3 emit/face/file air zhi data.beam beam_area_3 nevery 100 region circle3

################################################################################################

#dump                2 image all 10 image.*.ppm type type pdiam 0.2 surf proc 0.01 size 512 512 zoom 1.75 gline no 0.005
#dump_modify	     	2 pad 4

timestep            0.0001

stats		    	10
stats_style	    	step cpu np nattempt ncoll nscoll nscheck
run 		    	1000
Memory usage per proc in Mbytes:
  particles (ave,min,max) = 0 0 0
  grid      (ave,min,max) = 1.51379 1.51379 1.51379
  surf      (ave,min,max) = 0.00151062 0.00151062 0.00151062
  total     (ave,min,max) = 1.5153 1.5153 1.5153
Step CPU Np Natt Ncoll Nscoll Nscheck 
       0            0        0        0        0        0        0 
      10 0.00054311752        0        0        0        0        0 
      20  0.001023531        0        0        0        0        0 
      30 0.0015032291        0        0        0        0        0 
      40 0.0020711422        0        0        0        0        0 
      50 0.0026595592        0        0        0        0        0 
      60 0.0031607151        0        0        0        0        0 
      70 0.0036454201        0        0        0        0        0 
      80 0.0041637421        0        0        0        0        0 
      90 0.0046315193        0        0        0        0        0 
     100 0.0093636513     6286        0        0        0        0 
     110  0.012694359     6289        0        0        0       32 
     120  0.015004635     6289        0        0        0       24 
     130  0.019343615     6290        0        0        0       32 
     140  0.021728992     6289        0        0        0       24 
     150  0.026362181     6290        0        0        0    50064 
     160  0.044263363     6290        0        0        0    50320 
     170  0.065513134     6289        0        0        0    50304 
     180  0.081920147     6289        0        0        0    50304 
     190   0.10725307     6289        0        0     6175    99952 
     200   0.12774086    12643        0        0        0    52256 
     210   0.14250231    12644        0        0        0    51640 
     220   0.15965724    12645        0        0        0    50056 
     230   0.18949986    12646        0        0        0    41296 
     240   0.22037888    12641        0        0        0    28632 
     250   0.23025799    12607        0        0        0    69096 
     260   0.24839401    12418        0        0        0    63176 
     270   0.26464438    11925        0        0        0    59160 
     280    0.2824831    11162        0        0        0    56584 
     290   0.30549693    10305        0        0     6259   105352 
     300   0.32424331    15729        0        0        0    56296 
     310   0.33958936    14972        0        0        0    54240 
     320   0.35494518    14395        0        0        0    52056 
     330   0.37404156    13958        0        0        0    42736 
     340   0.38899612    13608        0        0        0    29312 
     350   0.39974165    13309        0        0        0    69360 
     360   0.41604209    12904        0        0        0    63000 
     370   0.43217325    12259        0        0        0    59160 
     380   0.45074725    11440        0        0        0    56712 
     390   0.47405624    10491        0        0     6209   104760 
     400    0.4948976    15784        0        0        0    55520 
     410   0.51079631    15023        0        0        0    54072 
     420   0.52559662    14419        0        0        0    51736 
     430   0.54373956    13903        0        0        0    42664 
     440    0.5585115    13537        0        0        0    29552 
     450    0.5696907    13223        0        0        0    68832 
     460   0.58570361    12835        0        0        0    62536 
     470    0.6051743    12227        0        0        0    58536 
     480   0.62315941    11384        0        0        0    55776 
     490   0.64428377    10461        0        0     6142   103680 
     500   0.66468477    15862        0        0        0    55368 
     510   0.67983079    15090        0        0        0    53648 
     520   0.69899869    14444        0        0        0    51384 
     530   0.72017717    13941        0        0        0    42520 
     540   0.73557281    13577        0        0        0    29816 
     550   0.74980855    13267        0        0        0    69496 
     560   0.76614285    12920        0        0        0    63144 
     570   0.78931165    12296        0        0        0    59216 
     580   0.80543399    11483        0        0        0    56664 
     590   0.82998824    10584        0        0     6234   105000 
     600   0.85163474    15998        0        0        0    55856 
     610   0.86665559    15201        0        0        0    54184 
     620   0.88438797    14523        0        0        0    51832 
     630   0.90215325    14029        0        0        0    42920 
     640   0.91428614    13669        0        0        0    29696 
     650   0.92806172    13342        0        0        0    69696 
     660   0.94757009    12955        0        0        0    62824 
     670   0.96777844    12367        0        0        0    58728 
     680   0.98766184    11536        0        0        0    56064 
     690    1.0160923    10604        0        0     6226   104496 
     700    1.0376849    15985        0        0        0    55824 
     710    1.0532875    15169        0        0        0    54032 
     720    1.0685585    14519        0        0        0    51736 
     730    1.0851488    14014        0        0        0    42560 
     740    1.0998807    13610        0        0        0    29312 
     750    1.1124499    13299        0        0        0    68952 
     760    1.1366379    12922        0        0        0    62736 
     770    1.1643436    12303        0        0        0    58832 
     780    1.1906183    11457        0        0        0    56152 
     790    1.2480645    10557        0        0     6218   104416 
     800    1.2747021    15907        0        0        0    55528 
     810    1.2901878    15082        0        0        0    53936 
     820    1.3106878    14426        0        0        0    51624 
     830    1.3297637    13957        0        0        0    42536 
     840    1.3444161    13597        0        0        0    29928 
     850    1.3573828    13308        0        0        0    69592 
     860    1.3802855    12946        0        0        0    63048 
     870    1.4042332    12317        0        0        0    58808 
     880    1.4257078    11493        0        0        0    56360 
     890    1.4526272    10593        0        0     6199   104528 
     900     1.478554    15904        0        0        0    55600 
     910    1.4942465    15040        0        0        0    54120 
     920    1.5092494    14383        0        0        0    51736 
     930    1.5260525    13871        0        0        0    43040 
     940    1.5402579    13464        0        0        0    29192 
     950    1.5503271    13162        0        0        0    68184 
     960    1.5689831    12758        0        0        0    62176 
     970    1.5851247    12155        0        0        0    58016 
     980      1.60165    11389        0        0        0    55552 
     990    1.6325927    10435        0        0     6074   102704 
    1000    1.6592443    15709        0        0        0    54672 
Loop time of 1.6593 on 4 procs for 1000 steps with 15709 particles

MPI task timing breakdown:
Section |  min time  |  avg time  |  max time  |%varavg| %total
---------------------------------------------------------------
Move    | 0.66323    | 0.84053    | 1.0813     |  18.3 | 50.66
Coll    | 0          | 0          | 0          |   0.0 |  0.00
Sort    | 0          | 0          | 0          |   0.0 |  0.00
Comm    | 0.047845   | 0.060282   | 0.072315   |   4.4 |  3.63
Modify  | 0.00040221 | 0.017757   | 0.037506   |  13.1 |  1.07
Output  | 0.0086641  | 0.0093712  | 0.011386   |   1.2 |  0.56
Other   |            | 0.7314     |            |       | 44.08

Particle moves    = 11032143 (11M)
Cells touched     = 11191026 (11.2M)
Particle comms    = 53795 (53.8K)
Boundary collides = 0 (0K)
Boundary exits    = 47302 (47.3K)
SurfColl checks   = 44025592 (44M)
SurfColl occurs   = 56650 (56.6K)
Surf reactions    = 0 (0K)
Collide attempts  = 0 (0K)
Collide occurs    = 0 (0K)
Reactions         = 0 (0K)
Particles stuck   = 0
Axisymm bad moves = 0

Particle-moves/CPUsec/proc: 1.66217e+06
Particle-moves/step: 11032.1
Cell-touches/particle/step: 1.0144
Particle comm iterations/step: 1.811
Particle fraction communicated: 0.00487621
Particle fraction colliding with boundary: 0
Particle fraction exiting boundary: 0.00428765
Surface-checks/particle/step: 3.99067
Surface-collisions/particle/step: 0.005135
Surf-reactions/particle/step: 0
Collision-attempts/particle/step: 0
Collisions/particle/step: 0
Reactions/particle/step: 0

Surface reaction tallies:
  id adsorb_test_ps2 style adsorb #-of-reactions 5
    reaction all: 124
    reaction C(b) --> C(g): 124

Particles: 3927.25 ave 5213 max 2629 min
Histogram: 2 0 0 0 0 0 0 0 0 2
Cells:      2 ave 2 max 2 min
Histogram: 4 0 0 0 0 0 0 0 0 0
GhostCell: 6 ave 6 max 6 min
Histogram: 4 0 0 0 0 0 0 0 0 0
EmptyCell: 6 ave 6 max 6 min
Histogram: 4 0 0 0 0 0 0 0 0 0
Surfs:    12 ave 12 max 12 min
Histogram: 4 0 0 0 0 0 0 0 0 0
GhostSurf: 0 ave 0 max 0 min
Histogram: 4 0 0 0 0 0 0 0 0 0
